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' We report the results of a search for an emission line from radiatively 

Q ' decaying dark matter in the ultra-faint dwarf spheroidal galaxy Willman 1 

^ ■ based on analysis of spectra extracted from XMM-Newton X-ray Observatory 

^ ■ data. The observation follows up our analysis of Chandra data of Willman 1 

I ' ( Loewenstein &: Kusenkd 120101 ) that resulted in line flux upper limits over the 

2 ■ Chandra bandpass and evidence of a 2.5 keV feature at a significance below the 

c/3 ■ 99% confidence threshold used to define the limits. The higher effective area of 

a ■ 

I— I' the XMM-Newton detectors, combined with application of recently developing 

,__! ■ methods for extended-source analysis, allow us to derive improved constraints on 

^ ■ the combination of mass and mixing angle of the sterile neutrino dark matter 

^ ' candidate. We do not confirm the Chandra evidence for a 2.5 keV emission line. 

(n: 

in . 

rn : 

O ' 1. Introduction 
(N ■ 

'^ . Sterile neutrinos represent a plausible dark matter candidate amenable to observation 

Sp ' ^^^ X-ray spectroscopy of locations in the universe where the dark matter surface density 

;-i . is high. Sterile neutrinos, gauge-singlet fermions that emerge in extensions of the Standard 

- - - Model of particle physics that explain ordinary neutrino masses, may comprise some or all 
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of the dark matter if the Majorana masses are below the electroweak scale (JKusenkd |2009| ) . 
If their mass lies in the 1-30 keV range, non-resonant oscillatior is produce sterile neutrinos 
at an abundance comparable to that inferred for dark matter (IDodelson fc Widrowl Il994[ ) 
- although several other production channels have been suggested (lKusenkdl2009l ). Sterile 
neutrinos in this mass range may explain the observed velocities of pulsars by anisotropic 
emission of sterile neutrinos from a cool ing neutron star born in a supernova explosion 
(iKusenko fc Segrg 119971 : Fuller et al.l 120031 ) and can facilitate early st ar formation (|Kusenko 



200 



9I). As a form of warm dark matter (WDM), sterile neutrinos ( Abazaiian et al. 



2001 



Petrakil 120081 : iBoyanovskyI l2008l : iBoyanovsky fc Wul 1201 ll : iDunstan et al.l l201l[ ) may resolve 
some of the discrepancie s (Ferrero et al. 2011, and references therein) betwe en Cold Dark 
Matter (CDM) models Jhovell et al.l boili iMenci. Fiore. fc Lamastral I2OI2I ) and the ob- 
served structure in the universe. A number of particle physics models have been pro- 
posed t o accommodate a sterile neutrino with the requisite mass and mixing, inc luding 
z/MSM ( Asaka fc Shaposhnikovl l2005l: iBovarskv. Ruchayskiv. fc Shaposhnikovi 12009), split 



seesaw (IKusenko. Takahashi. fc Yanagidall2010l ) , Higgs singlet (IKusenkd 



2006; 



Petraki fc Kusenko 



2008 ). and other models. For example, the split seesaw mechanism (IKusenko. Takahashi. fc Yanagida 

2OIOI ) with one extra dimension causes the mass and mixing of a sterile neutrino originating 

at a high scale to be exponentially suppressed in the low-energy effective theory, which makes 

it a good dark matter candidate, while preserving both the standard seesaw explanation for 

the neutrino ma sses and the explanation of t he matter-antimatter asymmetry of the universe 

by leptogenesis (JFukugita fc Yanagidalll986l ). 



The one-loop decay of relic keV sterile neutrinos into an active neutrino and photon 
with energy E^^ = mst(? /2 produces an X-ray emission line with width corresponding to the 
velocity dispersion of the dark matter distribution within the observed solid angle. While 
the prospe cts of discovering a dark - matter sterile neutrino in a labor atory present a daunting 
challenge (jAndo fc Kusenkd l2010l : iBezrukov fc Shaposhnikovl l2007l ) , the narrow decay liri e 
allows one to search for relic sterile neutrinos using X-ray telescopes (JAbazajian et al.ll200ll ). 
The high dark matter concentrations of dwarf spheroidal galaxies make them prime targets 
for dark matter searches, with the additional advantage of the total absence of competing 
intrinsic X-ray sources. Moreover, as the most dark matter dominated objects known they are 
not subject to systematic uncertainties that derive from decomposing the mass into baryonic 
and non-baryonic components - although assumptions about membership and dynamical 
equilibrium come into play in converting measured velocity dispersion profiles into dark 
matter mass estimates. We initiated a dedicated search for sterile neutrinos with the Suzaku 
X-ray telescope, placing new limits on sterile neutrinos over the 1-20 keV mass range from 
XIS s pectra of the Ursa Minor dwarf spheroidal galaxy (JLoewenstein. Kusenko. fc Biermann 
2009h . 
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We continued the search with Chandra observation of the ultra-faint dwarf spheroidal 
Willman 1 (Loewenstein & Kusenko 2010, ApJ, 700, 426; hereafter LKIO), considered at 
the time to be a particularly compelling target for dark matter searches (though its status 
subsequently changed somewhat; see below). 99% confidence line flux upper limits over the 
0.4-7 keV Chandra bandpass were derived and mapped to an allowed region in the sterile 
neutrino mass-mixing angle plane consistent with the Suzaku constraints. In addition we 
reported evidence for an emission line with flux below this threshold from radiative decay of 
a 5 keV sterile neutrino with a mixing angle in the narrow range where oscillations produce 
all of the dark matter and for which sterile neutrino emission from the cooling neutron stars 
can explain pulsar kicks (LKIO). This tentative result is best confirmed (or refuted) by 
utilizing the large effective area of the XMM-Newton EPIC detectors; and, we were awarded 
observing time on Willman 1 for this purpose. Here we present results of our analysis of 
these data, deriving more stringent constraints on sterile neutrino parameters and failing to 
confirm the best-fit estimate of the Chandra 2.5 keV line flux. 



2. Data Analysis 

2.1. Observation and Data Processing 

Willman 1 was observed with XMM-Newton in three segments in late October 2010: 
ObsID 0652810101 (for 29356 s on 10/22), ObsID 0652810301 (36054 on 10/25), and ObsID 
0652810401 (36218 s on 10/31). In addition, a region offset by 1 degree was observed (ObsID 
0652810201; 33856 s on 10/30). The last was primarily meant as a control in the event of 
a positive emission line detection; and, as none was found we will not consider these data 
further (although they were reduced and analyzed in parallel to the on-source data). We 
utilize data from the two EPIC-MOS (hereafter, EMOS = EMOSl + EM0S2), as well 
as the EPIC-PN (hereafter, EPN), CCD detectors. Data reduction, background modeling, 
and spectral extraction are conducted with the XMM-Newton Extended Source Analysis 
Software (XMM-ESAS) - part of the XMM-N ewton Science Ana l ysis S yste m (SAS vll.l is 



used t hroughout) - and methods as detailed in lKuntz fc SnowdenI (|2008[ ) and lSnowden et al. 



(120081 ). Although various extended source analysis techniques are adopted in X-ray dark 
matter searches, we are the first to apply this specfic methodology optimized for extended 
low surface brightness X-ray emission. A summary of our particular application of these 
techniques follows; further details on the approach and individual procedures and may be 
found in the XMM-ESAS "cookbook'fl 



^http://heasarc. gsfc.nasa.gov/docs/xmni/xmnihp_xininesas. html 
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In order to apply the latest calibration files and products, rather than use pipeline- 
processed data, the unprocessed event files are reduced to create filtered event files using the 
composite tools emchain and mos-filter [epchain and pn-filter) for EMOS (EPN) datasets. 
This results in the removal of time intervals of prominent soft proton (SP) flaring - identified 
by hard (2.5-12 keV) energy band background count rates deviating by more than 1.5a from 
the mean of a Gaussian distribution that characterizes the count rate histogram. Out-of-time 
(OOT) EPN events that correspond to photons registered during CCD readout are processed 
in parallel to ordinary EPN events for subsequent correction (see below). Exposure times 
for the final cleaned event files are shown in Table 1. Primarily due to fiaring, data from 
~one-third of the EMOS, and ~ 60% of the EPN, exposures are excluded in the end. 



2.2. The Extraction of Source and Background Spectra 

We exclude background point source as follows. Source detection is conducted with 
the SAS vll composite edetect_chain tool as a standalone procedure and as incorporated in 
the XMM-ESAS cheese procedure. These are used in the construction of region masks in 
detector coordinates for use in spectral extraction. We considered masks that (a) minimize 
the excluded regions by restricting only the cores (the inner 8" in radius) of the brightest 
sources (those with likelihood> 40), (b) that minimize contamination from point sources 
(32" radii for sources detected with likelihood> 10), as well as (c) an intermediate choice 
(16" radii for sources detected with likelihood> 10). We also considered three spectral 
extraction regions. The first corresponds to 1.7x the half-light radius (1.7 x 25 pc = 232" 
at the 38 kpc Willman 1 distance) within which the mass is particularly well-determined 



(jAmorisco fc Evansll201ll ) and for which we apply mask "a" above, the second to 150 pc 



(818"; mask "b"), and th e third to 100 pc ( 545"; mask "c") - the mass distribution is 



estimated within 100 pc in lStrigari et al.l (|2008[ ). The latter two are of the same order as the 



optical extent of Willman 1, and may be compared to the tidal radius of 930 pc (Section 3) 
and a virial radius on the 10^ M© mass-scale of ~ 25 kpc. 

Table 1. Final Good Time Intervals in ks 



Obsid 


EMOSl 


EMOS2 


EPN 


0652810101 
301 
401 


15.37 
21.86 
27.30 


19.21 
23.13 
28.43 


9.53 

15.47 

16.16 


total 


64.53 


70.70 


41.16 
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Files were processed with FLAG == and screened to retain only those events with 
PATTERN <= 12 (four-or-fewer pixel events) for EMOS, and PATTERN <= 4 (single 
and double pixel events) for EPN, spectra as a means to exclude non- X-ray events. The 
XMM-ESAS mos-spectra and mos-back [pn-spectra and pn-back) procedures are used to ex- 
tract EMOS (EPN) source and quiescent particle background (QPB) spectra, and to compute 
spectral response files. This method utilizes filterwheel-closed data, data from the unexposed 
corners of archived XMM-Newton observations, and ROSAT All-Sky Survey (RASS) data 
to produce a model QPB spectrum. Chips found to be in anomalous states with elevated 
low-energy background according to the XMM-ESAS criteria (EM0Sl-ccd4 and -ccd5 for 
ObsID 0652810301 and ObsID 0652810401, EMOS2-ccd5 for ObsID 0652810101) are ex- 
cluded, and EM0S-lccd6 was no longer operating at the time of these observations. The 
detector extraction areas are calculated using the proton_scale task. The spectral extraction 
procedure generates instrument response matrices in the form of separate redistribution ma- 
trix (rmf) and effective area function (arf) files, taking into account the extended nature of 
the emission. 



2.3. Spectral Analysis 

2.3.1. The Baseline Spectral Model 

The XMM-ESAS procedures create model QPB background spectra appropriate to the 
conditions during the observation under consideration; however, a number of internal and 
astrophysical background components remain that must be included in any spectral model. 
This approach is both more accurate and more conservative than one where blank-sky 
background is subtracted, given the directional dependence of the intensity and spectral 
shape of the c osmic contribution, and the time dependence of the particle and instrumenta l 



backgrounds ( JKuntz fc SnowdenI l2008l : ISnowden et al.l l2008l : iLoewenstein &: Kusenkd l2010l ) . 



Moreover, in the event that there is dark matter line emission the "blank-sky" includes a 
contribution to the signal that ought not be subtracted. 

The instrumental background includes a number of fluorescent lines. The strongest of 
these - the 1.49 keV Al Ka (EMOS and EPN) and 1.75 keV Si Ka (EMOS only) features are 
sufficiently variable that the XMM-ESAS model QPB spectrum includes a smooth bridge 
over the relevant spectral region, and the lines must be included in spectral modeling in 
the form of narrow Gaussians. Other problematic EPN lines are avoided by restricting the 
bandpass to < 7 keV. 

Even after filtering, contamination from residual SP flaring and solar wind charge ex- 
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change (SWCX) are commonly present. These are modeled using a broken power-law and 
a pair of unresolved emission lines at 0.56 keV (OVIl) and 0.65 keV (OVIII) that dominate 
SWCX spectra, respectively. Since energy from the SP component is directly deposited in 
the detectors, diagonal response files are used for this component in place of the telescope 
responses described above. 

The astrophysical X-ray background includes multiple components that we model with 
thermal plasma (apec) models to account for the emission from the Local Hot Bubble 
(LHB) and the Milky Way Halo (MWH), and a power-law for the unresolved point sources 
mostly originating from background AGN (CXB). The e lemental compositioii in the apec 
components are set at their solar abundances as defined in lAsplund et al.l (|2009[ ). We use the 
HEASARC background tooOto extract the ROSAT All-Sky Survey (RASS) spectrum from a 1° 
radius aperture centered on the position of Willman 1 for use in constraining the parameters 
of these components. 



2.3.2. Fitting Techniques and Best-Fit Baseline Models 

Spectra are fitted using XsPEC version 12.7lj. The final total spectral model we adopt 
may be expressed, in XsPEC notation, as bknpower -|- gaussian -|- gaussian -|- con- 
stant X constant x (gaussian -|- gaussian -|- apec -|- (gaussian -|- apec -|- apec 
-|- constant x powerlaw) x TBabs), where bknpower represents the residual SP com- 
ponent (with its distinct diagonal response matrix), the first pair of gaussians the Al Ka 
and Si Ka instrumental lines, and the second pair of gaussians the SWCX emission. The 
remaining astrophysical background consists of the unabsorbed apec LHB and absorbed 
(two-temperature) apec -|- apec MWH, as well as the absorbed powerlaw CXB emission. 
The constant factor m ultiplying the CXB r epresents the departure of the CXB intensity 
from its all-sky average (JKushino et al.ll2002l ) due to cosmic variance and the fact that the 
resolved portion of the CXB is removed. The Tuebingen-Boulder ISM absorption model 
(jWilms. Allen, fc McCrayl |2000[ ) . tb abs, is applied with the column density fixed at the 



Galactic value of 1.13 x 10^° cm ^ (JDickey fc LockmanI Il990l ). The two multipliers that 



act on all of the sky components, constant x constant, represent the size of the spectral 
extraction region (see above) so that all normalizations correspond to fluxes per solid angle, 
and a factor to account for any detector calibration offset. 



^http://heasarc. gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg. pi 
■^ http://heasarc.gsfc. nasa.gov/docs/xanadu/xspec/ 
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There are datasets from three detectors for each of the three on-source ObsIDs. We 
initially analyze each dataset separately. As is usually the case, we find that the EMOSl 
and EM0S2 spectra are very similar for each observation, and we always analyze these 
in tandem with all LHB, MWH, and CXB component parameters tied together. Due to 
systematic differences between the EMOS and EPN CCDs, particularly with regard to details 
in particle background for extended sources, we always fit these data separately and derive 
independent constraints on sterile neutrino line emission from each set of detectors. Joint 
fits would invite the introduction of a bias in the best fit model parameters and could lead to 
the underestimation of errors by artificially degrading the model goodness-of-fit (e.g., Baldi 
et al. 2012). 

Although the corresponding systematics are less severe we choose not to coadd the 
datasets from the three separate observation intervals, fitting them first separately and then 
simultaneously as we shortly describe in detail. Our final results are taken from the latter. 
All parameters are tested to see if they are well-determined, if fits are sensitive to their 
values, and if limits on additional line emission over a range of representative energies are 
significantly affected by their variation. When thus justified, these are fixed as described 
below. The following is based on extensive experimentation, from which models with a 
minimum number of significant, variable parameters emerge that lead to a robust limit on 
sterile neutrino line emission. 

In the simultaneous fits, the Al Ka and Si Ka energies and widths are fixed at the 
values determined from the individual observation fits. The SWCX line energies are fixed at 
0.56 and 0.65 keV, and their widths at 1 eV (negligible compared to the EPIC resolution). 
The LHB temperature is fixed at 0.15 keV, the MWH temperatures at 0.10 and 0.34 keV 



based on fits that included the RASS spectrum, and the CXB slope at 1.46 ( jKushino et al. 



2002|). The EM0S1/EM0S2 offset is set to unity. 



In what follows, we restrict discussion to the 100 pc radius circular aperture spectral 
extraction region. The increase in counts provided superior statistical accuracy to the 42.5 
pc (1.7 X the half-light radius) region. Given the systematic mass uncertainties (see below), 
we judge that this compensates for the formally smaller uncertainty in the mass enclosed 
within the smaller aperture. While the larger extraction aperture for the 150 pc region 
encompasses a larger flux, the effect of vignetting increases the ratio of non-X-ray to X-ray 
events and degrades the statistical accuracy. 

We first consider the unbinned, unsubtracteco spectra over the "global" 0.3-10 (0.4-4) 



^Here, "background-subtracted" and "unsubtracted" always refer to whether or not the QPB - only - 
is subtracted prior to fitting, and "unbinned" refers to the default energy intervals defined by the detector 



keV energy range for the EMOS (EPN), fit by minimizing the modification of the C-statistic 



( ICashl I1979I ) implemented in XsPEC as cstat (we refer to this approach as ^^cstat-nbs^^ , 
where nbs refers to "no background subtraction"). Since the CXB is much smaller than 
the QPB, its normalization (relative to the all-sky average; see above) is also fixed at values 
determined from the separate observation fits. The SWCX norms are set to 0, and the relative 
observation-to-observation offsets are set to 1. Thus the free parameters are the separate SP 
broken power-law slopes for each detector and each observation, and the SP broken power- 
law energy break for each detector (linked at an identical value across observations). The free 
normalizations include those for the SP component and instrumental lines for each detector 
and each observation, and those for the LHB and MWH (tied for each set of EMOS or EPN 
detectors for each observation). 

The following minor statistical issue arises with respect to the EPN regarding out-of-time 
(OOT) events. The XMM-ESAS pn-back task produces an OOT-subtracted spectrum that 
is not purely Poissonian, and the cstat statistic that we sometimes employ is strictly valid 
only for data with Poisson errors. We address this by including an explicit OOT background 
in the se cases, scaling the OOT spectra cre ated by the XMM-ESAS pn-spectra task by 



0.063 (JNevalainen. Markevitch. fc Lumbll2005l ). Results using the presubtracted spectra are 



consistent, as expected given the small departure from Poissonian, with the approach we 
adopt - although the latter yields somewhat better fits. 

Although these are not utilized to derive best-fits or sterile neutrino line flux constraints, 
for purposes of illustrating the relative strengths of the spectral model components and the 
overall signal-to-noise of the data we show the coadded EMOS (EMOSl and EM0S2, all 
observations) and EPN (all observations) total spectra and mean best-fitting model in Figure 
1. 

Because of the prominence of the (unsubtracted) instrumental fluorescent lines, deriving 
limits on additional line emission is highly problematic in the 1.2-1.9 (1.3-1.65) keV energy 
range based on EMOS (EPN) spectra. As a result it is sensible to split the spectrum into 
low- and high-energy segments on either side, and treat these separately. Moreover, since 
some of the model components are negligible in one or the other of these segments one can 
apply simpler models with parameters that are a subset of those in the global fit. Ultimately, 
the constraints on sterile neutrino line emission are derived from fits to these low- and high- 
energy sub-spectra, with the role of global fits purely to provide a means of setting these 
up. 

We proceed with the analysis described above in the low energy spectral segments (0.3- 



channel boundaries. 



-9- 




1 2 

Energy (keV) 




Energy (keV) 



Fig. 1. — The left panel (a) {right panel (b)) shows the total co-added EMOS (EPN) unsubtracted spectrum 
(black errorbars), and mean best-fitting model (red curve) composed of astrophysical (sum of LHB, MWH, 
and CXB; blue curve) background, particle background (SP; green broken power-law), and instrumental 
fluorescent lines (orange). This is meant to illustrate overall S/N for the unbinned spectra (note that the 
EMOS energy channels are wider than the EPN channels) and the spectral decomposition; co-added spectra 
are not used to derive spectral analysis results. 



1.2 keV for the EMOS, 0.4-1.3 keV for the EPN) by freezing the instrumental hne norms 
at the values derived from the global fits, reactivating the SWCX lines (with normalizations 
tied across observations), and with the SP component model converted to a single-slope 
power-law allowed to vary from observation to observation. We follow with fitting of the 
high energy (1.9-10 keV for the EMOS, 1.65-4 keV for the EPN) spectral segments by freezing 
the LHB, MWH, and SWCX normalizations and converting the SP component back to a 
broken power law for the EMOS; parameters (slopes, break energies if applicable) are tied 
across observations. Spectra and best-fit models are shown in Figure 2. 

In our second approach we consider (QPB) background-subtracted, binned (minimum 
15 counts per bin) spectra fit by minimizing the x^ statistic (we refer to this approach as 
"X^-^s")- The global EPN bandpass may now be extended to 7 keV (residual instrumental 
features at ~ 4.5 and ~ 5.5 keV that preventing this in the unsubtracted spectrum are now 
subtracted out). With the increased prominence of the CXB relative to the SP, we now fix the 
SP parameters (though not their normalizations) to the best-fit values determined above, 
and thaw the CXB norm (relative to its all-sky value) which is tied across observations. 
For the low-energy sub-spectra, the instrumental line normalizations are again fixed at the 
values derived from the global fits and the SWCX lines are reactivated. In addition, the 
CXB normalization is now fixed. In following this with fitting of the high energy sub-spectra 
the LHB, MWH, and SWCX normalizations are frozen, but the CXB normalization thawed 
and allowed to separately vary for each observation. Again, for illustrative purposes, the 
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Fig. 2. — The left panel (a) {right panel (b)) shows the 6 (3) individual EMOS (EPN) spectra (with 
some binning) and best-fit models for the simultaneous cstat fits to the unsubtracted spectra. Low- and 
high-energy segments, although shown together here, are separately fitted (see text for details). 



coadded subtracted spectra and corresponding mean best-fitting model are shown in Figure 
3. Spectra and best-fit models for the individual spectra in each of the two energy intervals 
are shown in Figure 4. 




1 2 

Energy (keV) 




1 2 

Energy (keV) 



Fig. 3. — Same as Figure 1 for the subtracted spectra (with some binning). 

Finally, we also consider joint fitting of the total (unsubtracted) and background (QPB) 
spectrum. Here the background spectrum, which is not Poissonian, is binned to a minimum 
of 50 counts per bin and fitted by application of x^ statistics, while the total spectrum is 
unbinned and fit by applying the cstat statistic. Thus, the best-fit model parameters are 
determined from minimization of a hybrid statistic. The global fits proceed as in the first 
approach with SP and instrumental parameters and normalizations tied in the models for 
the total and QPB spectra, and the other (SWCX, LHB, MWH, and CXB) component 
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Energy (keV) Energy (keV) 

Fig. 4. — Same as Figure 2 for the subtracted spectra (with some binning, and on a logarithmic scale). 

normalizations set to zero in the latter (thus no additional parameters are introduced). For 
the low-energy sub-spectrum, the SP is characterized by a broken power-law fixed across 
observations and detectors. This is also the case for the high-energy EMOS spectrum, while 
the SP is assumed to be a single-slope power-law for the high-energy spectral EPN segment. 

The fit statistics for all three approaches are displayed in Table 2, along with the 90% 
(A — statistic = 2.71, where "statistic" refers to the fit statistic - either x^, the modified 
C-statistic, or the hybrid statistic described above) confidence upper limits on the surface 
brightness of a narrow emission line with energy fixed at 2.5 keV. Negative emission line 
fiuxes are permitted when deriving these limits. 



2.4. Emission Line Flux Limits at 2.5 keV and Other Energies 

Having established the baseline models as described above, we derive upper limits as in 
Loewenstein et al. (2009), and LKIO. An unresolved Gaussian component, stepped in 10 eV 
intervals over the relevant bandpass, is added to the baseline model, and A — statistic = 9.21 
(99%) upper confidence levels on the line fiux (that is permitted to be negative) are computed. 
Limits for the low- and high-energy segments are separately computed with the spectral 
model parameters fixed or variable as described above - except for the SWCX component 
that we always fix at their best-fit values (with the - necessary - inclusion of this component 
in XMM-Newton spectra, limits obtained in the ~0. 55-0. 66 keV region using this method are 
provisional). The limits, in 150 eV bins for the x^-bs and cstat-nbs approaches applied both 
to the EMOS and EPN detectors, with zooms on the higher-energies for the latter, are shown 
in Figure 5. Given the similarity of the limits from the two approaches we henceforth adopt 
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Table 2. Fit Statistics - 100 pc Extraction Region 



method 


bandpass 


detector 


Stat 


best-fit-stat 


Sx(2.5) 


cstat-nbs 


global 


EMOS 
EPN 


cstat 


3993/3841 
2108/2153 






lo-en 


EMOS 
EPN 




340/337 
492/532 






hi-en 


EMOS 




3342/3225 








EPN 




1416/1403 


1 1/1 + 1-51 


X^-bs 


global 


EMOS 
EPN 


x' 


1545/2488 
2032/2025 






lo-en 


EMOS 
EPN 




221/341 
499/500 






hi-en 


EMOS 




1021/1866 


0-70+l'25 






EPN 




1354/1328 


1 70+1.48 
l-'>^-1.41 


hybrid 


global 


EMOS 
EPN 


hstat 


4359/4535 
2812/2505 






lo-en 


EMOS 
EPN 




390/442 
654/681 






hi-en 


EMOS 




3630/3819 


n qq+o-so 






EPN 




1732/1608 


9 Yo+1.46 



Note. — Shown arc the fit statistics (per degree-of-frccdom) for the 
best-fit model determined using the three analysis approaches - unbinned, 
background-unsubtracted source spectral fitting with the (modified) C- 
statistic C^cstat-nbs" ) , binned, background-subtracted source spectral fitting 
with the x^ statistic ("x^-^''^"); joint unsubtracted-source/background spec- 
tral fitting using a combination of the x^ s-nd C-statistic ("hybrid") - see 
text for details. Also shown are the best-fit and 90% confidence limits on the 
average surface brightness in units of 10^* photons cm^^ 
narrow emission line with energy fixed at 2.5 keV. 
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the x^-bs limits. The x^-bs constraints are formally more restrictive for the EPN detector in 
the 4-7 keV bandpass. However as this is a result of an overall shift in the allowed line flux 
and not of an improvement in accuracy and, given the residual background artifacts in this 
energy range (§2.3.2), we also restrict our subsequent discussion to constraints derived from 
the EMOS detectors. 



The expected line flux from sterile neutrino radiative decay in Willman 1, produced by 
active-sterile neutrino transitions with standard assumptions about the lepton asymmetry 
and thermal history of the universe below ~ 1 GeV (LKIO and references therein; also, see 
below), is shown by the solid lines in Figures 5ab (the star in (c) shows the prediction for 
rrist = 2.5 keV). The broken line shows the allowed range taking the impact of hadronic 
uncertainties, which affect the relation between the mass and the mixing angle, into account 
(jAsaka. Laine. fc Shaposhnikovl 120071 ) . An average total (sterile neutrino) dark matter sur- 
face density of 200 M© pc~^ is adopted (see below). A sterile neutrino produced in this way 
in sufficient abundance to compose all of the dark matter produces an emission line that ex- 
ceeds the observed limits if ?7ist > 5 keV. That is, the presence of a 2.5 keV line corresponding 
to the parameters of the best-fit Chandra estimate is not confirmed. We further illustrate 
this in Figure 6 that compares, for both approaches and for the EPN and EMOS, measured 
spectra and best-fit models that include a narrow 2.5 keV emission line with strength fixed at 
that predicted for radiative decay of sterile neutrinos produced by non-resonant oscillations. 
There is no evidence of such a feature in the XMM-Newton spectra. 
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Fig. 5. — Upper limits on the emission line surface brightness, in units of photons cm~^ s^^arcmin^ 
and averaged over the inner 100 pc of Willman 1, for both EMOS and EPN detectors and x^'bs and 
cstat-nbs approaches as indicated in the embedded legend {left panel (a)), with focus on the high-energy 
spectra segment {x^-bs approach) in middle (b) and right (c) panels. The solid and broken lines show the 
expected line flux from sterile neutrino radiative decay in Willman 1, for an average total (sterile neutrino) 
dark matter surface density of 200 Mq pc^^ (see text), taking into account the hadronic uncertainties in 
production (Asaka et al. 2007). The star in (c) shows the mean prediction for nist = 2.5 keV. 
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Limits on Sterile Neutrino Parameters from the XMM-Newton Spectrum of 

Willman 1 



The equations relating the dark matter projected surface mass density Sdm, X-ray ob- 
servables (line energy E^, and line surface brightness Eune), and sterile neutrino parameters 
(mass r/ist, mixing angle 9, and fraction of dark matter in sterile neutrinos /st), are the 
following (see, e.g., LKIO and references therein): 



m. 



st 



2E. 



11 



r.. 



^IVa 



5.52 X 10" 



-32 



sin^^ 



10 



-10 



VkeVy 



and 



-^line 



3.95 X lO^^r.^^^^J, 



-yUa J st 



-'dm 



Mopc- 



E^ 



-1 



photons cm" 



'^ s ^arcmin ^, 



(2) 



(3) 



where the expression for the decay rate of relic keV sterile neutrinos into an active neutrino 
and photon, T^^^^y^, is given for Majorana sterile neutrinos. 

As in LKIO we estimate the ave rage dark matter surface mass density withi n 100 pc 

based on projecting the best-fit NFW (JNavarro. Frenk. fc Whitdll997l ) mass model inlStrigari et al 



(|2008[) : however, we now introduce a (tidal) truncation radius of 930 pc (ISanchez-Conde et al. 
2011 ) by adopting the "n = 2 BMO" generalization of the NFW profile ( jOguri &: Hamana 
201 ll ). This yields 135 M© pc~^ from the dark matter in Willman 1, to which we (conserva- 



tively) add 65 Mq pc~^ (LKIO) associated with dark matter in the Milky Way halo to obtain 
a fiducial total line-of-sight surface mass density SHmflOO Pc) = 200 MfD pc~^. The resultin g 
mass profile is consistent with other estimates ( IWolf et al.l l2009l : lAmorisco fc Evansll201ll ). 
However Sdm(100 pc) is uncertain by a factor of 2 or more, considering the errors associated 
with those in the velocity dispersion profile, as well as the dynamical state of Willman 1 
(see discussion below). The form of the mass distribution is also uncertain and, as would 
be expected for warm dark matter, there is evidence that the dark matter density profile 
is flatter and smoother than might be expected in a CDM universe on the relevant mass 
scale (Gilmore et al. 2007; Cholis & Salucci 2012 and references therein) and may even be 
constant over 100 pc. The impact of these uncertainties on the average density is smaller 
than that of these other sources (unlike the case for dark matter annihilation calculations 
where the emissivity depends on the square, rather than the first power, of the density). 

The upper limit on the dark matter radiative decay rate as a function of energy from 
equation (3) with /st = 1 is shown in Figure 7, where we also show the limits for Sdm(100 pc) = 
100 M0 pc~^ (or, equivalently, /st = 0.5). As in LKIO and Loewenstein et al. (2009) we map 
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this into excluded regions in the rrist-O sterile neutrino parameter space under each of the 
following two separate assumptions: (1) that all of the dark matter is composed of sterile 
neutrinos produc ed by some unspecified mec hanism, (2) that non-resonant oscillations, as 



first suggested by iDodelson fc Widrowl (J199J), produced sterile neutrinos at the abundance 



determined by mst and 9 as calculated in Asaka et al. (2007) (Figure 8). The first region is 
where Sune, calculated from equations (l)-(3) assuming fst = 1, exceeds our inferred upper 
limits; the second where the flux from sterile neutrinos produced at the minimal abundance 
from oscillations exceeds these limits. The latter provides an absolute constraint, assuming 
only a standard early thermal history of the universe, since the oscillations cannot be turned 
off. 



3.1. Discussion and Conclusions 

We have derived upper limits on the radiative decay of dark matter into keV photons 
in the ultra-faint dwarf spheroidal galaxy Willman 1, and used these to place very general 
constraints on the mass and mixing angle of any sterile neutrinos that contribute to the dark 
matter. In doing so we adopted a total surface mass density (including Milky Way dark 
matter in the line of sight) of Sdm(100 pc) = 200 Mq pc~^ in the solid angle subtended 
by the XMM-Newton aperture corresponding to 100 pc at the Willman 1 distance of 38 
kpc. While this estimate is based on equilibrium dynamical models fit to the stellar velocity 
dispersion profile, the dynamical status - and even the very nature of Willman 1 - were re- 
examined subsequent to our original X-ray observation of Willman 1 with Chandra. Recent 
optical spectroscopy supports the case that Willman 1 is (or was) a dwarf spheroidal and not 
a star cluster; however - in addition to the limited statistics of the velocity dispersion profile - 
doubts as to whether Willman 1 is in dynamical equilibrium introduce a level of uncertainty in 
the Willm an 1 dark matter mass that must now be estimated from non-equilibrium dynamical 



modeling ( iWillman et al.ll201l[ ). Thus while our constraints are formally comparable to 
previous limits based on dwarf spheroidal X-ray spectroscopy (Loewenstein et al. 2009, 
Riemer-S0rensen & Hansen 2009, LKIO, and references therein), we cannot claim to provide 
additional new constraints on sterile neutrino parameters. 

On the other hand we applied (for the first time) recently developed extended-source 
analysis methods to demonstrate the robustness and sensitivity of the XMM-Newton EMOS 
and EPN detectors to weak line emission that may arise from dark matter radiative decay 
- despite the fact that ~half the data were discarded due to the effect of flares. Clearly 
XMM-Newton has an important role to play in such studies, with its larger effective area 
complementing the lower and more stable background of Suzaku and the superior spatial 
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resolution (crucial for sources that, unlike dwarf splieroidals, have prominent intrinsic discrete 
source emission) of Chandra. 

We find no confirmation of the Chandra evidence for an emission line at 2.5 keV 
(Figures 5-6). Moreover, the abundance of mgt = 5 keV sterile neutrinos produced by 
non-resonant oscillati ons has recently bee n strongly proscribed based on Chandra imaging 
spectroscop y of M31 (Watson et al.ll201ll ). albeit in a region where baryonic matter is still 
prominent. IWatson et al.l ( 201ll ) derive an upper limit m^t < 2.2 keV for sterile neutrinos 
produced by this mechanism at the abundance required to explain all of the dark matter 
(fst = 1). Combined with the lower bo und from phase-space considerations of ~ 0.3-0.4 keV 
( Tremaine &: GunnllQTQt Angus 2010 ), as we ll as those at higher energy from the Milky Way 
and CXB ( lYiiksel. Beacom. fc Watson! |2008| ) . this implies that accounting for dark matter 
in this manner with sterile neutrinos may remain viable only at masses where the decay line 
lies in a spectral region dominated by blended emission features from the LHB and MWH - 
and, thus, sensitivity to detection at X-ray CCD spectral resolution is limited. As a result, 
the most significant expansion of the searchable parameter-space boundary requires, above 
all, the sort of leap in spectral resolution of diffuse sources possible with microcalorime- 
ter arrays such as the Soft X-ray Spectrometer that is part of the Astro-H observatorjO 
scheduled for launch in 2013. As investigation of this lower rUst regime proceeds, the region 
where off-res onant oscillations produce pulsar kicks at the high end of the observed velocity 
distribution ( JFuller et al.l 120031 : iKusenko. Mandal. fc Mukherjed 120081 ) will be probed. 



Alternatively, the sterile neutrino abundance may be subdominant {fst < 1), or ac- 
counted for by means in addition to the DW production mechanism. The production of sterile 
neutrin os via oscillations rnay be resonantly enhanced in the p r esence of nonzero lepton asyrn - 
metry (jShi fc FuUerl Il999l : iKishimoto. Fuller, fc SmithI 120061 : iLaine fc Shaposhnikovl |2008| ) , 
thus opening up a larger region in the mgt — 6 plane consistent with observed line emission 
constraints and /st = 1- Other suggested production mechanisms, i.e. from inflaton or Higgs 
decay, do not involve oscillations and are therefore independent of the mixing angle ( IKusenko 
20061 : IShaposhnikov fc Tkachevll2006t iPetraki &: Kusenkd 120081). Furthermore, if sterile neu- 



trinos are produced above the electrow eak scale (JKusenkd l2006t iPetraki &: Kusenkd 12008 



Kusenko. Takahashi. &: Yanagidall2010f) , or if some other sterile neutrinos decay and produce 
entropy ( lAsaka et al.l l2006t Fuller et al.ll2009l . l201l[ ). both the abundance and the clustering 
properties of dark matter are affected. However should the pulsar kick and non-resonant os- 
cillation domains be ruled out, the expected detectability of sterile neutrino radiative decay 
lines becomes essentially undetermined. 



^ http://heasarc.gsfc. nasa.gov/docs/astroh/ 
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Despite ongoing dedicated production and direct detection experiments, as well as 
intensive indirect searches, successful discovery of CDM in the form of weakly interact- 
ing massive particles (WIMPs) remains elusive. Given the prevalence of sterile neutrinos 
in extensions of the standard model that explain the generation of neutrino masses, the 
eme rgence of the keV scale as a natural on e for some re cently proposed production scenar- 
ios (JKusenko. Takahashi. fc Yanagidall2010l : lMerlell2012l ). and indications that cosmogonies 
dominated by WDM may explain disc repancies between CDM and observations of small 



scale structure an d galaxy formation (JLovell et al.l l201ll : iMenci. Fiore. fc Lamastral 12012 



Yue fc Chenll2012l ). strong physical motivation for their search persists. In addition, there 
are other decaying keV dark matter candidates, one of which - m oduli dark matter - we are 
currently investigating ( JLoewenstein. Kusenko. &: Yanagidall2012[ ). 
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Fig. 6. — Unbinned EMOS {upper left panel (a)) and EPN {upper right panel (b)), and binned EMOS 
{lower left panel (c)) and EPN {lower right panel (d)) spectra (black errorbars) with best-fit models (and 
their decomposition) that include a narrow 2.5 keV emission line with strength fixed at that predicted for 
radiative decay of sterile neutrinos produced by non-resonant oscillations (see Figure 5). The upper and 
lower curves correspond to the cstat-nbs and y^-bs approaches, respectively. The color coding follows Figure 
1, with the additional line component in orange. 
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Fig. 7. — Excluded dark matter radiative decay rates (filled region) as a function of line energy derived 
from the upper limit on the average emission line surface brightness emerging from the inner 100 pc of 
Willman 1, assuming I]dm(100 pc) = 200 Mq pc~^ of decaying dark matter. The yellow histogram shows 
the corresponding limits for Edm(100 pc) — 100 M0 pc~^ - most of which would originate in the Milky Way 
halo. 
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Fig. 8. — The shaded region in the rUst — sterile neutrino parameter plane is generally excluded assuming 
only the standard cosmological history below the temperature where production by neutrino oscillations 
occurs. The region to the right of the solid line is excluded if all of the dark matter is composed of sterile 
neutrinos produced by some (unspecified) mechanism. 



